Unstable spectrum of a relativistic electron beam 
interacting with a quantum collisional plasma: 
application to the Fast Ignition Scenario 

A. Bret, F. J. Man'n Fernandez, J.M. Anfray 

ETSI Industriales, Universidad de Castilla-La Mancha, 13071 Ciudad Real, Spain 
E-mail: antoineclaude . bretSuclm . es 



Unstable spectrum of a relativistic electron beam. 



2 



Abstract. Quantum and coUisional effects on the unstable spectrum of a relativistic 
electron beam-plasma system are investigated through a two-fluids model. Application 
is made to the near target center interaction of the relativistic electron beam in the 
Fast Ignition Scenario. Partial degeneracy effects are found negligible while the most 
influential factors arc the beam temperature and the electron-ion collision frequency 
of the plasma. The introduction of the latter triggers some oblique unstable modes of 
much larger wave length than the coUisionless ones. Transition from the coUisionless 
regime to the resistive one is thus documented and found discontinuous. 



PACS numbers: 52.57.Kk,52.35.Qz 
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1. Introduction 

The Fast Ignition Scenario (FIS) for Inertial Confinement Fusion assumes that a 
pre-compressed Deuterium- Tritium target is ignited by a Petawatt Laser generated 
relativistic electron beam, originated near the border of the pellet [H E]- This 
scenario offers a very interesting theoretical setting to plasma physics, as it includes 
the interaction of a relativistic electron beam with a plasma which density varies over 
several orders of magnitude. When the electron beam reaches the center of the pre- 
compressed pellet, the plasma density is so high that quantum effects are to be accounted 
for. Interestingly, if beam-plasma instabilities are not welcome at the beginning of 
the beam travel to the core, they may not be so deleterious at the end of it because 
they contribute to the beam stopping and energy deposition. Much efforts have been 
devoted so far to the beam-plasma interaction physics involved in the whole process 
[3 m El E], but the investigation of quantum effects is still in its infancy. As far as 
collisional effects are concerned, recent works focused on the filamentation instability 
(see Refs. [TJ [8] and references therein). While this instabihty is expected to govern 
the system at the beginning of the beam path because of the high beam to plasma 
density ratio [9], the low density ratio near the target center would rather point towards 
an "oblique modes" driven system, where the dominant mode wave vector makes an 
oblique angle with the beam flow. Regarding these modes, quantum effects have never 
been assessed, collisional ones are almost unknown and the interplay between these two 
remains unidentified. Finally, it has been now established that beam temperature is a 
key parameter [3], HOl |TT] when describing the unstable spectrum. 

Our intent is to construct a theory accounting for all the aforementioned effects 
in order to identity the most influential ones, and focus on them in later works. The 
model presented here includes therefore the following features: 1) a relativistic beam 
with transverse and parallel temperatures, 2) a weakly degenerate and 3) collisional 
plasma, and finally, 4) the calculation of the whole unstable spectrum in order to spot 
the most unstable modes in this setting. We thus implement a quantum collisional 
two-fluids model, as previous works [12] demonstrated that such formalism is able to 
reproduce the results from a waterbag kinetic theory [TOl [11]. Note that on the one 
hand, the fluid treatment of thermal effects requires sub-relativistic temperatures (see 
remark after Eq. [2]) while on the other hand, thermonuclear burn of the target should 
deeply modify the physics considered in the sequel. The process considered here is thus 
the beam interaction with the pre-ignited and pre-compressed FIS target core. 

2. Basic equations 

We consider the interaction of an homogenous and infinite electron beam of density 
and velocity v^q with a plasma of density Upo- A plasma return current at velocity 
VpQ cancels exactly the beam current with nf,QVbQ = UpoVpQ. Ions are considered a fixed 
neutralizing background. Furthermore, our focus on the near target center interaction 
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implies n^Q <^ UpQ. Both beam and plasma electrons share the same conservation 
equation, 

^ + V ■ (n, V,) = 0, (1) 

but we use quite different Euler equations for the two species. 

On the one hand, we assume a relativistic beam which density is low enough for 
quantum effects to be neglected. We also neglect here all sort of collisions due to the 
relativistic velocity involved; collisionality is here restricted to the plasma species [TSllS]. 
The beam Euler equation reads therefore 

(v6 • Vjpb = -g E + , (2) 



dt \ c J Uf, 

where g > is the electron charge and p;, = 'ymvb the relativistic momentum. The 
kinetic pressure term will be dealt with according to the guidelines set in Ref . [12] : The 
pressure gradient is first expressed as VPb = SkBThVuf, [m [5] before a temperature 
tensor is introduced in the linearized equation (see Eqs. [HllQl and more comments bellow). 
Such an adiabatic treatment requires sub-relativistic temperatures [15l[l6]. While this 
is not a strong requirement for the plasma, such condition is more stringent for the 
beam. 

On the other hand, plasma electrons are drifting with the non-relativistic velocity 
Vpo = {ni^o / npQ)vu) <^ c. They suffer collisions with the background plasma and are 
weakly degenerate. Thus, we use for them the following modified Euler equation 

+ Vp ■ V Vp = - ^ E + ^ , - . vp 

at m \ c ' 



The first two terms are just the Lorentz and the friction force where v is the collision 
frequency between the plasma return current and the background ions. The third 
and fourth terms have been proposed recently as the quantum correction to the fluid 
equations for a finite temperature Fermi plasma [T^. They result from the combination 
of the so-called Bohm pressure term [TH [191 120] for a completely degenerate jellium, 
with a correction proportional to the thermal velocity v'^p = kBTp/m accounting for the 
finite plasma temperature Tp. Equation ([3]) can be seen as a combination of Eq. (19) 
in Ref. [E] and Eq. (19) in Ref. [20] where spin and relativistic effects are neglected. 
The function G in the equation above reads [T7] . 

where fi is the chemical potential, = kBTp and Li^ is the polylogarithm function, 

CO ^ 



=1 
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Denoting Ep the Fermi energy of the electron gas, G ~ ^(3Ep for fc^^p ^ Ep. For 
Bp < 2kBTp, 

1 /Ep^'/' 



G ~ 1 + —= . (6) 

A core density of 10^^ cm~^ yields Ep = 0.78 keV. With a core temperature of a few 
hundreds of eV, we find that neither the classical ksTp ^ Ep nor the cold jellium 
approximations kBTp -C Ep are really valid for the electron gas because its temperature 
is quite close to its Fermi energy. Partially degenerate thermal corrections to the Bohm 
pressure are thus required in Eq. ([3]). 

The beam is relativistic but classical while the return current is non-relativistic but 
quantum. We can thus account for both relativistic and quantum effects through Eqs. 
f l2|3|) and do not need quantum relativistic corrections in Eq. ([3]). 

In order to compute the dispersion equation, Eqs. ffT]|2][51) are linearized assuming 
small perturbations of the equilibrium variables proportional to exp(ik ■ r — iujt) where 
i"^ = —1. In the present relativistic diluted beam regime, we expect oblique modes to 
govern the system [9]. We thus consider both parallel and perpendicular components of 
the wave vector to make sure the fastest growing modes are not overlooked. Choosing 
the z axis as the beam direction, we can set k = {kx, 0, kz) without loss of generality. 
The linearized equations are 

k ■ Vol /„\ 

for the perturbed beam and plasma densities uji, and for the beam 



im7b(k • Vfoo - t^) 



7,2 



El 



c 



- 3^— A;sT, ■ k, (8) 



ribo 

where T;, is the temperature tensor 

( Tb± \ 
Tb= Tb^ . (9) 
V Tfeii J 

Although some covariant [211 ESI EHl El] or kinetic [HI ES] treatment of the beam 
temperature would be more appropriate because of the relativistic regime involved here, 
it was proved in Ref. [12] that the present formalism can reproduce the results from 
a waterbag kinetic theory, providing the perpendicular beam temperature parameter 
Tb± is re-scaled by a factor 1/^/%. At any rate, both waterbag kinetic theory and 
the present model are quite limited in their treatment of temperature. Nevertheless, 
waterbag theories are quite common j26| [3l [10] as a first approach to kinetic effects, 
especially in the relativistic regime. The linearized Euler equation for the plasma reads, 

im(k ■ Vpo - uj)\pi = - q ^Ei + XeL^^^ _ mu(ypo + Vpi) (10) 

4m^ Upo ^Hpo 
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From this stage, Eqs. (ITIISIITU]) are solved to express the perturbed velocities v^i and 
Vfci in terms of Ei and Bi. The first order magnetic field is then eliminated through 
Bi = (c/co')k X El, yielding an expression of the current, 

J = ^ njoVjo + rijoVji + rijiVjo, (11) 

j=b,P 

in terms of Ei only. Finally, the equation above is inserted into a combination of 
Maxwell-Faraday and Maxwell- Ampere equations, 

4?7r 

— k X (k X El) + El + Ji(Ei) = 0, (12) 

UJ"^ UJ 

which gives the dielectric tensor. This tensor has been here symbolically computed with 
the Mathematica Notebook described in Ref. [27j and is provided as "Supplementary 
Material" to the paper. It is expressed in terms of the following dimensionless variables, 
nto ry kt;6o n ^2 ^kBTb±,\\ 

where = AitUp^q^ / m is the plasma electronic frequency, and 

where the parameter G is given by Eq. (j4]). While Gc measures the strength of the 
zero temperature Bohm pressure term, Bt accounts for finite temperature effects. In 
this respect, the bridge between the completely degenerate and the classical cases is 
established through this parameter [T7] . 



3. Quantum and thermal effects 

We start from the most simple cold-beam/cold-collisionless-plasma case for which the 
maximum growth rate 5m is reached for an oblique wave Zm vector with [28] , 

— = 7^ {—] ^zM ~ 1, Z^M > 1- (15) 

In the limit of this cold regime, there is not one single most unstable mode but a 
continuum of oblique modes growing at the same rate [29j. Because we deal with 
multiple effects at once, we focus on some typical beam and plasma parameters for FIS, 
and look at the influence of beam temperature, plasma degeneracy and collisionlality 
on the most unstable mode as well as its growth rate. We thus consider from now on 
parameters values extracted at peak compression from a global simulation performed 
by Ren et. al. [30] , 

rife = 10^2 cm-^ 7b = 3, T^xji = 100 keV, 

Up = 10^^ cm-^ Tp = 1 keV, u = OAup. (16) 

Other simulations consider the same values for the core parameters [36]. The fixed 
dimensionless parameters in Eqs. (JT3ll are thus, 

a = 10-^ pb±,\\ = 0.76, Gc = 1.3 x 10"^ Bt = 6.5 x 10"^ (17) 
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(A) (C) 




Figure 1. (Color online) Growth rate normalized to the cold coUisionless one given by 
Eq. (|15ll8p . in terms of Z = 'kvbo/i^p- (a) Cold classical coUisionless case, (b) Former 
case plus Bohm pressure term, (c) Former case plus quantum thermal corrections, 
(d) Thermal effects included, plus all quantum corrections. Collision frequency is zero 
(r = 0) for all figures. 

while the maximum growth rate from Eq. (|T5l) is 

Sm = 0.022a;p. (18) 

3.1. Quantum effects 

Starting from the cold, classical coUisionless growth rate map on Fig. W^a), we 
successively add "cold" quantum effects on Fig. [I](b), then partial degeneracy correction 
on Fig. [U^c) and beam temperature on Fig. [Hd). It is clear from Fig. [^b) that the 
Bohm pressure term hardly changes the global picture whereas Fig. [D^c) shows more 
influence from the partial degeneracy correction, although the maximum growth rate is 
not affected. 
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Quantum thermal corrections arising from the "G" term in Eq. ([3]) can be compared 
with the so-called Bohm pressure term. Considering the unstable spectrum emphasizes 
wave vectors with k ~ Up/ c, and comparing the two quantum terms, the Bohm term is 
found secondary while 

Qc < Qt, (19) 

In view of the numerical values given by Eqs. ( JT71) . the condition above is largely fulfilled. 
The plasma is therefore hot enough for the Bohm pressure term to be negligible. The 
two quantum parameters would be comparable for Gt about 5 x 10^ smaller. Because 
6t oc oc Tp, this would require Tp ~ 1/50 eV. 

The fluid quantum theory of the filamentation instability [311 [32j showed that 
quantum effects reduce small wave length instabilities through quantum interferences. 
We find here a similar trend all over the unstable spectrum. By lowering the growth 
rates at large kx (i.e., k±), quantum effects alone are found to single out one most 
unstable mode where the classical cold case yields a continuum of them. The most 
unstable mode on Fig. [l](c) is 

Z^M ~ 1-06, ~ 4.07, with — = 0.021 (20) 

UJp 

which growth rate is close to the one of the cold classical system given by Eqs. fllSflSp . 
Although quantum corrections are found to have a noticeable effect on the unstable 
spectrum, they do not alter the fastest growing mode. 



3.2. Beam temperature effects 

Beam temperature effects are very important, as evidenced on Fig. [l](d). In accordance 
with some recent conclusions drawn from a more elaborated kinetic relativistic model [5], 
the fastest growing coUisionless mode is here oblique. Also, most of the unstable modes, 
filamentation included [33l [3ll |3], have been stabilized and the unstable spectrum is 
now restricted to a very narrow region of the Z plane. The maximum growth rate has 
been reduced to 15% of its cold classical value, and the fastest growing mode is now, 

Z^M ~ 1.35, ZxM ~ 1.08, with — = 0.0033. (21) 

Up 

Note that the same calculation canceling all quantum effects gives the very same result. 
So far, plasma degeneracy is found to affect poorly the dominant unstable mode. 

The growth rate map evolves dramatically when collisions are accounted for. As 
observed on Fig. [2|^a-c), collisions strongly reduce the former largest growth rate. But 
the most remarkable evolution is the arising of oblique coUisional unstable modes at 
much lower k. Indeed, the fastest growing mode and its growth rate is now, 

Z^M ~ 0.0078, ZxM ~ 0.029, with — = 4.7 x 10"^ (22) 

Up 

CoUisional effects are therefore essential in the present setting, and we now focus on 
them. 
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Figure 2. (Color online) Growth rate of the whole unstable spectrum in terms of 
Z = hvb/ujp including all effects and for various collision frequencies. Note the 3 
orders of magnitude drop of the maximum growth with respect to Fig. [T]^a) . 



4. Collisional effects 

The arising of collisional modes more unstable than the collisionless ones has already 
been mentioned in the literature for the filamentation instability [5j. Still for the 
filamentation instability, the shift of the fastest growing mode to lower fc's has also 
been reported [H [H [35] . Indeed, such a shift is desirable if the transition is to be made 
from the collisionless regime, where the filaments transverse size is of the order of the 
plasma skin depth, to the resistive regime llHl [36], where filaments transverse size is 
about the much larger beam skin depth (at least in the diluted beam regime). What is 
evidenced here is the extension of this trend to the oblique modes which indeed govern 
the diluted relativistic regime [9]. 

We can monitor between Figs. [H^d) and [2|(a-c) the "birth" of the collisional mode 
at much lower k, which eventually surpasses the collisionless fastest growing mode when 
T is increased. On the one hand, collisions strongly limit the growth of the collisionless 
oblique modes while on the other hand, triggering some low k's unstable modes. 

We now focus on the transition dynamic. How does the system switches from 
one dominant mode to the other? Figs. [2](a,b) display the unstable spectrum for 
intermediate collision frequencies between Fig. [It^d), where r = 0, and Fig. [2](c), 
with r = 0.4. Comparison suggests the transition happens here for r < 0.1. A finer 
evaluation yields r = 0.085, which is much smaller than the 0.4 value considered here for 
FIS. Therefore, and according to the present model, the end of the beam path definitely 
pertains to the collisional regime, in the sense that unstable modes stemming from 
coUisionality govern the system. Noteworthily, the transition is abrupt. We do not have 
one mode evolving continuously from one regime to the other. Instead, as collisions 
are progressively "switched-on" , dominant collisional unstable modes are progressively 
stabilized in one region of the unstable spectrum, while in a completely different region, 
collisional modes grow faster and faster. When some critical collision frequency is 
reached (here u = 0.085cc;p), the latters overcome the formers and the dominant mode 
"jumps" from the resonant part of the spectrum ~ 1 to a much lower one. Figures 
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Table 1. Most unstable wave vector and its growth rate (in cOp units) accounting for 
the various effects presented. Except in the first case, there is always one single most 
unstable mode. 



Model 


3m 


ZzM 


ZxM 


Cold only 


0.022 


1 


> 5 


Bohm term 


0.022 


1 


13 


All Quantum terms 


0.021 


1.06 


4.07 


Quantum + thermal 


0.0033 


1.35 


1.08 


Thermal only 


0.0033 


1.35 


1.08 


All effects + t = 0.4 


4.7 X 10-'' 


0.0078 


0.029 



[2]^b,c) show that the new modes are shifted down by some 2 orders of magnitude along 
the Zz and axis. Remembering the beam to plasma density ratio is here 10~^, two 
orders of magnitude also separate the plasma skin depth from the beam one. The low 
k modes found here generate patterns of the "correct" typical size, namely, the beam 
skin depth one. 

Table [H allows for a global picture of the results gathered in Sees. [3] & HI Every 
effect added to the cold case calculation reduces the maximum growth rate. Quantum 
effects only slightly do so, while the beam thermal spread dramatically slows down 
the fastest growing mode. This most unstable mode remains on the very same branch 
on the dispersion equation as long as collisions are neglected. By setting r 7^ 0, a 
new branch appears at much lower k, which overcomes the collisionless branch from 
T = 0.085. Furthermore, the hot and collisional system is about 3 orders of magnitude 
less unstable than its cold counterpart. 

5. Conclusion 

By implementing a relativistic fluid model for beam plasma interaction, we could assess 
the relative importance of collisions and quantum effects for the instabilities arising 
when a hot relativistic electron beam interacts with the dense core of a pre-ignited 
and pre-compressed FIS target. We consider the beam as relativistic, collisionless 
and classical. The plasma return current is degenerate, non-relativistic and collisional. 
Also, we computed the whole unstable spectrum in order to spot the fastest growing 
unstable mode, whether it is two-stream, filamentation or oblique-like. For the typical 
FIS parameters mentioned in Eqs. (|T6l) . Table [1] summarizes the growth rate of the 
fastest growing mode as well as its location accounting for the various effects. 

Quantum effects result from the inclusion of the zero temperature Bohm pressure 
term, and of a quantum thermal correction. The Bohm term hardly modifies the 
unstable spectrum and the largest growth rate, as can be checked from Table [1] or Figs. 
[Tl^a,b). However, by damping large k± unstable modes, it singles out one most unstable 
wave vector. This trend amplifies when including quantum thermal corrections. The 
normal component of the most unstable wave vector decreases even more, but the value 
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of 6m remains unchanged. Finally, calculations accounting for beam temperature and 
collisions are insensitive to the inclusion of any quantum effects. Although some kinetic 
treatment could be needed to settle the case, it seems that thermal and collision effects 
are definitely the dominant ones in FIS context whereas plasma degeneracy can be 
neglected. Furthermore, plasma temperature should rise as the beam starts depositing 
its energy, resulting in a even less degenerate plasma. 

Including beam temperature brings about a major modification of the unstable 
spectrum (see Fig. [Hi). As expected from previous kinetic studies |3l|3l], filamentation 
modes are here completely damped. Even oblique modes close to filamentation have 
been stabilized and by virtue of the very small beam to plasma density ratio, the 
unstable spectrum has been reduced to a thin line starting from ~ 1. Both the 
location and the growth rate of the most unstable mode are affected, and the typical 
size of the generated patterns is the plasma skin depth. 

Finally, and this may be the most important conclusion of this paper, the inclusion 
of a strong collision electron-ion frequency for the plasma return current triggers some 
very low k oblique unstable modes. These longer wavelength modes overcome the 
collisionless ones as soon a.s u > Uc = 0.085ujp. Note that such modes also exist 
along the perpendicular direction, but the dominant ones are oblique. Given the high 
collision frequency considered in our case (0.4ti;p), the collisional regime should definitely 
be switched-on. The transition from one regime to the other is discontinuous when 
crossing the Uc threshold because the two competing modes are located in remote places 
of the spectrum. A discontinuous evolution of the growth rate was already pointed out 
at z/ = by Molvig, in a kinetic analysis of the magnetized filamentation instability 
|33j . In the present case, the growth rate is a continuous function of k and u, and 
the discontinuity arises only when considering the maximum growth rate over the k 
spectrum. Additionally, the discontinuity occurs here for = z^c 7^ 0. 

By growing at much smaller k (see Tabled]), the collisional mode generate patterns 
which typical size is now the much larger beam skin depth, reminiscent of the so-called 
resistive filamentation instability [I3l[36]. At the present stage, we think more studies are 
needed to identity the present collisional modes with this instability, but similarities are 
striking. A previous cold model [37] accounting for the same collision frequency for both 
the beam and plasma electrons, failed to generate such modes. Furthermore, we checked 
that the present model yields fast growing collisional filamentation modes in the absence 
of beam temperature. This latter factor thus seems responsible for the obliqueness of 
the fastest growing mode, while the plasma collisionality (together with the beam non- 
collisionality) seem to be the reasons for their very existence. It is probable that a 
certain degree of collision is acceptable in the beam before the collisional regime vanishes. 
Since the expression of the dielectric tensor is known, it should be possible to access 
the collisional modes analytically, gaining thus a deeper understanding of the transition 
uncovered here. Appropriate approximations are currently being developed allowing to 
simplify the dielectric tensor. Finally, beam temperature effects on the collisional modes 
need to be kinetically assessed to go beyond the present sub-relativistic temperature 
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limit. 
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